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We study the use of truncated normal-ordered three-nucleon interactions in ab initio nuclear structure cal- 
culations starting from chiral two- plus three-nucleon Hamiltonians evolved consistently with the similarity 
renormalization group (SRG). We present three key steps: (i) a rigorous benchmark of the normal-ordering 
approximation in the importance-truncated no-core shell model (IT-NCSM) for ■*He, "'O, and '*''Ca; (ii) a direct 
comparison of the IT-NCSM results with coupled-cluster calculations at the singles and doubles level (CCSD) 
for '*0; and (iii) first applications of SRG-evolved chiral NN-f3N Hamiltonians in CCSD for the medium-mass 
nuclei '^^'•o and '*'''***Ca. We show that the normal-ordered two-body approximation works very well beyond 
the lightest isotopes and opens a path for ab initio studies of medium-mass and heavy nuclei with chiral two- 
plus three-nucleon interactions. At the same time we highlight the predictive power of chiral Hamiltonians. 
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Two decades of experience with ab initio nuclear struc- 
ture calculations, including the pioneering work with the 
Green's Function Monte Carlo (GFMC) approach HI, have 
shown that three-nucleon (3N) interactions play an impor- 
tant role for understanding the structure of nuclei systemat- 
ically from first principles. Recent advances on nuclear in- 
teractions derived within chiral effective field theory (EFT) 
put additional emphasis on the consistent inclusion of real- 
istic 3N interactions in nuclear structure calculations i^-IH. 
However, for most large-scale many^body approaches, such 
as the no-core shell model (NCSM) ^ or the coupled-cluster 
(CC) method 101, the full inclusion of 3N interactions in- 
creases the computational cost by orders of magnitude and 
often renders calculations impossible that are routinely per- 
formed with two-nucleon (NN) Hamiltonians. In order to re- 
solve this dilemma — the need for 3N interactions versus the 
tremendous computational cost — one might resort to approx- 
imate schemes for including the 3N interaction, particularly 
when aiming at a controlled approximation rather than an ex- 
act solution of the many-body problem. Along these lines, 
density-dependent NN interactions are being used for nuclear 
matter and finite nuclei to simulate the effects of 3N inter- 
actions in a computationally simple approximation (for re- 
cent applications see Refs. Isl-llOll). A related but more rigor- 
ous and improvable scheme to derive an approximate lower- 
rank form of the 3N interaction is based on a normal-ordered 
Hamiltonian with respect to a nucleus-specific reference state. 

In this letter we present rigorous benchmarks and sys- 
tematic applications of the normal-ordering approximation 
for chiral 3N interactions. We use the importance-truncated 
NCSM (IT-NCSM) iH [H for the ab initio solution of the 
many-body problem for '^He, '^O, and ^"Ca and compare 
ground-state energies and expectation values obtained with 
the full and the truncated normal-ordered 3N interaction. We 
then apply the normal-ordered 3N interaction in CC calcula- 
tions for the ground-states of medium-mass nuclei up to '***Ca 
and discuss the implications for nuclear structure predictions 
with chiral NN-I-3N Hamiltonians. 



Normal-Ordering Approximation. Transforming the 
many-body Hamiltonian into a normal-ordered form with 
respect to a nontrivial vacuum state is a standard technique in 
quantum many -body physics. Using creation and annihilation 
operators, a'y and a,,, for a single-particle basis |v) defined 
with respect to a trivial zero-body vacuum state, the operator 
of a 3N interaction formally reads 
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where V^^J^^J^^ - {yiV2V3\V2N\l^if^2f^3) are antisymmetrized 
matrix elements, and the operator Aj^'J/, ' = fli,flvT " ■ 
is a short-hand for a normal-ordered product of creation and 
annihilation operators with respect to the trivial vacuum (vac- 
uum normal-ordering). 

Instead of the trivial vacuum, we can choose a reference 
state |<l>i-ef) being a Slater-determinant of A single-particle 
states and reinterpret the creation and annihilation operators 
as particle or hole creation and annihilation operators with re- 
spect to this new vacuum. We have to rearrange the creation 
and annihilation operators in the particle-hole picture to estab- 
lish normal-ordering with respect to |<l>ref) (reference normal- 
ordering). Using Wick's theorem and A^JJ^^ ■ as a short-hand 
for the reference normal-ordered product, we obtain for the 
3N interaction 
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for the zero- 
for the one-body (IB) 
for the two-body (2B) term, and 



T,V|ff2a3 



KlJf^A = KIm2m\ for *e residual three-body (3B) term, where 
a, labels the occupied single-particle states in |Oief). 
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The interesting aspect of the reference state as compared to 
the trivial vacuum is that it contains information about the spe- 
cific many-body system under consideration — the reference 
state provides a first approximation to, e.g., the ground-state 
of a closed-shell nucleus. Based on this information, contri- 
butions of the 3N interaction are demoted to lower-particle 
ranks. This is the basis of the normal-ordered n-body (NOnB) 
approximation. 

Benchmark of the Normal-Ordering Approximation. To 
quantify how well such a truncation works, we perform ab 
initio calculations for the ground states of closed-shell nu- 
clei in the IT-NCSM. The IT-NCSM allows us to include the 
exact 3N interaction just as well as any NO«B approxima- 
tion. The underlying Hamiltonian contains the chiral NN in- 
teraction at N^LO of Ref. 



3|] and the local chiral 3N inter- 
action at N^LO of Ref. |ll4t|. The low-energy constants co 
and ce are taken from a fit to the ground-state energy and 
the yS-decay half-life of A = 3 systems 11511 . This initial 
Hamiltonian is transformed through a similarity renormaliza- 
tion group (SRG) evolution at the two- and three-body level 
to enhance the convergence behavior of the many-body cal- 
culation lll6l4l8ll . The SRG evolution represents a continu- 
ous unitary transformation parametrized by a flow-parameter 
or, with a - Q corresponding to the initial Hamiltonian. We 
will mainly consider two types of SRG-evolved Hamiltoni- 
ans: The NN+3N-full Hamiltonian starts with the initial chi- 
ral NN-I-3N Hamiltonian and retains all terms up to the 3N 
level in the SRG-evolution, the NN+3N-induced Hamiltonian 
omits the chiral 3N interaction from the initial Hamiltonian, 
but keeps all induced 3N terms throughout the evolution. The 
3N terms in both Hamiltonians have quite different character- 
istics, which makes them useful for benchmarking the normal- 
ordering approximation. In addition we will employ a range 
of values of the flow parameter a to generate an even larger 
set of test cases. 

We start with a direct comparison of IT-NCSM calculations 
for the ground-state energies of ^He and '^O using either the 
exact 3N interaction or the N02B approximation. The ref- 
erence state is always the QhQ. Slater-determinant composed 
of harmonic-oscillator single-particle states. The IT-NCSM 
energies as function of the model-space truncation parameter 
Nmax are presented in Fig.[T] The comparison of the converged 
values of the ground-state energy obtained with the N02B ap- 
proximation and with the exact 3N interaction reveals a mul- 
tifaceted picture. The largest relative deviation at the level of 
2% is observed for '^He with the NN+3N-induced interaction. 
The N02B approximation leads to a overbinding of about 0.4 
MeV for the softest Hamiltonian {a - 0.08 fm"*) and 0.6 MeV 
for the hardest Hamiltonian {a - 0.04 fm**). This systematics 
is completely reversed when using the NN-i-3N-full Hamil- 
tonian, here the softest interaction produces the largest devi- 
ation of about 0.3 MeV, whereas for the hardest interaction 
the N02B approximation coincides with the exact 3N cal- 
culations. For '^O we observe deviations below 1% for the 
ground-state energy, again with a nontrivial dependence on 
the Hamiltonian. For NN+3N-induced the N02B approxima- 
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FIG. 1: (color online) IT-NCSM ground-state energies for ''He and 
"'O as function of A'n,i„ for the NN+3N-induced and the NN+3N- 
full Hamiltonians for a range of flow parameters: a = 0.04 fm"* (•), 
0.05 fm* (♦), 0.0625 fm* (a), 0.08 fm^* (■). Solid symbols con-e- 
spond to the exact 3N interaction, open symbols to the N02B ap- 
proximation. Error bars indicate the uncertainties of the threshold 
extrapolations of the IT-NCSM. Data points are connected by straight 
lines to guide the eye, beyond the largest an exponential extrap- 
olation fitted to the last four data points is shown. 



tion consistently overestimates the binding energy by about 1 
MeV, for NN+3N-full the N02B approximation overbinds by 
less than 1 MeV for a - 0.08 fm"* and underbinds by less than 
1 MeV for a = 0.04 fm'*. 

For a comprehensive picture of its anatomy, we analyze 
the expectation values of the 3N interaction at different lev- 
els of the NO«B approximation using IT-NCSM eigenstates 
obtained with the exact 3N interaction for "*He, '''O, and "^"Ca 
for fixed A^max- Figure |2] summarizes these expectation val- 
ues of the 3N interaction for a set of NN+3N-induced and 
NN+3N-full Hamiltonians. For '^O and '^'^Ca a similar pattern 
emerges: The N02B approximation does reproduce the ex- 
pectation value of the exact 3N interaction very well, both for 
the NN+3N-induced and the NN+3N-full Hamiltonian. The 
pattern observed for the sequence of NO«B approximations is 
different for both types of Hamiltonians. ForNN+3N-induced 
the IB and 2B contributions of the normal-ordered Hamilto- 
nian have opposite sign, with the IB contribution being sig- 
nificantly larger, whereas for the NN+3N-full Hamiltonian the 
IB and 2B contributions are both attractive and of similar size. 
In all cases the OB contribution is the largest and overestimates 
the exact 3N expectation value. For "^He the pattern is differ- 
ent. The OB term does not provide the largest contribution and 
underestimates the 3N expectation value. The signs and rela- 
tive sizes of the IB and 2B terms again depend on the Hamil- 
tonian and the N02B approximation still shows a sizable de- 
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FIG. 2: (color online) Anatomy of the NOnB approximation of the 
ground-state energies of ''He, '^O, and '*''Ca. The bar charts show the 
expectation values of the 3N interaction computed at different levels 
of the normal-ordering approximation, i.e., NOOB, NOIB, N02B, 
and exact 3N. We employ the NN+3N-induced and NN+3N-full 
Hamiltonians, each with two values of a (see labels). We use the 
eigenstates obtained for the exact 3N interaction in = 10 for 
-•He and "^O and iV,^ax = 8 for ""Ca, all at HQ. = 20 MeV. 

viation from the exact 3N expectation value — consistent with 
the discussion of Fig.[T] 

This case study shows that there is no universal pattern and 
no hierarchy in the individual NOnB contributions. The size 
of the individual terms and also the deviation of the N02B ap- 
proximation from the exact 3N result depends on the Hamil- 
tonian, the nucleus, and the oscillator frequency. Nonethe- 
less, the 3N expectation values in Fig.|2]and the ground-state 
energies in Fig. [1] demonstrate that the N02B approximation 
works very well beyond the lightest nuclei. 

Application in Coupled-Cluster Theory. After validating 
the N02B approximation, we are now applying it in ground- 
state calculations for heavier closed-shell nuclei in the frame- 
work of the coupled-cluster method. Coupled-cluster theory 
is a natural framework, since normal-ordering of the Hamilto- 
nian with respect to a reference state is inherent to the for- 
mulation of the approach. We have developed an efficient 
coupled-cluster code using the /-coupled scheme discussed 
in Ref. d, which enables us to go to very large model 
spaces. We limit ourselves to singles and doubles excita- 
tions (CCSD), which has been shown to be a good approxi- 
mation for soft SRG-evolved interactions lfl9ll . An additional 
approximation present in the CCSD calculations for techni- 
cal reasons is a truncation of the 3N matrix elements entering 
the N02B to harmonic-oscillator principal quantum numbers 
e\ +62+^3 < £3 max = 14. Unlike the A^^ax truncation of the 
IT-NCSM, the e^ax truncation of CCSD does not imply this 
3B energy cut so that we miss some 3N matrix elements for 

^max ^ 6- 

In a first step, we confront the CCSD results for '^O with 
the previous IT-NCSM results, both using the N02B approx- 
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FIG. 3: (color online) Comparison of the ground-state energies of 
'*0 obtained in IT-NCSM and CCSD including 3N interactions at 
the N02B level for the NN-l-3N-induced and the NN-l-3N-full Hamil- 
tonians with a = 0.04 fm* (•), 0.05 fm" (♦), 0.0625 fm"* (a), and 
0.08 fm'' (■). 

imation. Figure [3] shows the convergence of the ground- 
state energies in both methods using the NN-i-3N-induced and 
NN-i-3N-full Hamiltonian. We observe a very good agreement 
of the converged ground-state energies, with the IT-NCSM 
giving 1 to 2 MeV more binding. This difference is consis- 
tent with the contributions expected from triples coiTections 
and the missing 3N matrix elements with £3 max > 14. The 
latter point has been confirmed by comparing to lower £3 max 
cuts. Altogether, the CCSD calculations for '^O with soft 
SRG-evolved NN-I-3N Hamiltonians in N02B approximation 
provide a ground-state energy within 1 to 2% of the IT-NCSM 
results with the exact 3N interaction. 

Using CCSD with the N02B approximation we can now 
study the systematics of ground-state energies with SRG- 
evolved chiral NN-I-3N Hamiltonians beyond '^O. Follow- 
ing the analysis of Ref. Iil61 we discuss the a-dependence 
observed with the NN-only, the NN-i-3N-induced, and the 
NN-i-3N-full Hamiltonians for '^O and ^*0, shown in Fig. 
m and for "^''Ca and "^'^Ca, shown in Fig. |5] For all nuclei 
we observe the same pattern: The NN-only Hamiltonian ex- 
hibits strong a-dependence of the converged ground-state en- 
ergies hinting at induced 3N interactions. Their inclusion at 
the NN-i-3N-induced level eliminates the a-dependence, thus 
demonstrating that induced 4N contributions originating from 
the initial NN are irrelevant. The converged energies, there- 
fore, correspond to the solutions for the initial chiral NN in- 
teraction. We obtain -120.2 (-t-0. 8) MeV for ""O ground- 
state energy, -152.1 (-1-0.5) MeV for -^Q, -343 {+6) MeV for 
4"Ca, and -392 {+!) MeV for ^^Ca using the NN-i-3N-induced 
Hamiltonian at a = 0.04 fm"* for Cmax — 14. The numbers 
in parenthesis give the change when going to a = 0.08 fm"* 
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FIG. 4: (color online) CCSD ground-state energies for and -''O 
as function of fimax for the three types of Hamiltonians (see column 
headings) using the N02B approximation for a range of flow param- 
eters: a = 0.04 fm* (•), 0.05 fm"* (♦), 0.0625 fm"* (a), and 0.08 fm-* 
(■). The filled symbols for the NN+3N-full Hamiltonian are for the 
standard chiral 3N interaction with cutoff 500 MeV, the open sym- 
bols for a modified 3N interaction with cutoff 400 MeV (see text). 
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In addition to the standard chiral 3N interaction f]3\ with 
cutoff momentum of 500 MeV, we also employ a chiral 3N 
interaction with a modified cutoff of 400 MeV and C£ = 0.098 
refitted to reproduce the "^He binding energy. We keep the 
value cd = -0.2 as in the standard 3N interaction. Effec- 
tively the lower cutoff reduces the strength of the two-pion 
terms of the 3N interaction and limits them to lower momenta. 
As a result the a-dependence and thus the induced 4N con- 
tributions are reduced significantly. This allows for a quan- 
titative comparison of the NN-i-3N-full predictions with ex- 
perimental binding energies. We obtain ground-state energies 
of -126.4 (-1.9) MeV for "'O, -164.8 (-2.8) MeV for ^^Q, 
-357 (-6) MeV for ^OCa, and -403 (-8) MeV for ^'^Ca using 
a = 0.04 fm^ with the change when going to c = 0.08 fm^ 
given in parenthesis. The agreement with experiment is re- 
markable. For and -"^O the predictions based on the chi- 
ral NN-t-3N Hamiltonian are consistent with the experimental 
binding energies. Even for "^"Ca and ^^Ca the agreement with 
experiment is surprisingly good, given the fact that no infor- 
mation beyond "^He was used to fix the Hamiltonian. This 
is evidence that this chiral NN-I-3N Hamiltonian, although 
not fully consistent regarding the chiral order and the cutoff 
choice in NN and 3N terms, contains all the relevant physics 
for nuclear structure predictions over a large mass range. 

Conclusions. We have demonstrated that the N02B ap- 
proximation allows for accurate nuclear structure calculations 
using SRG-evolved chiral NN-I-3N Hamiltonians in cases, 
where the inclusion of the exact 3N interaction is computa- 
tionally too demanding. Therefore, it provides a valuable tool 
to exploit the full physics potential of chiral Hamiltonians. In 
this context we have shown that a chiral 3N interaction with 
reduced cutoff can yield binding energy systematics consistent 
with experiment — ongoing investigations of the spectroscopy 
of p- and sd-shell nuclei confirm the quality and universality 
of this Hamiltonian. A generalization of the normal-ordering 
approximation to open-shell systems and multi-determinantal 
reference states is subject of present research. 
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FIG. 5: (color online) Same as Fig.|4]for *''Ca and "'^Ca. 



as a measure for the residual a-dependence. These results 
are in very good agreement with the CC results reported in 
for the bare chiral NN interaction. 
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When including the initial 3N interaction, i.e., when using 
the NN-i-3N-full Hamiltonian, the ff-dependence reemerges, 
indicating that 4N terms induced by the initial 3N interaction 
become sizable. These CCSD results confirm the findings of 
Ref. d and extend the systematics to heavier nuclei. 
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